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Abstract 

The first observation of the B j? —> D*^K± decay is reported using 3.0fb _1 of proton- 
proton collision data collected by the LHCb experiment. The D* T mesons are 
reconstructed through the decay chain D U —> 'yDj (K^K^ir^). The branching 
fraction relative to that for B j? — > D*~n + decays is measured to be 

B(B° S ^ D*^K ± )/B(B° S ^ D*~n + ) = 0.068 ± 0.005^;^, 

where the first uncertainty is statistical and the second is systematic. Using a recent 
measurement of B(B® — > D*~n + ), the absolute branching fraction of B® —> D*^K± 
is measured as 

B(B® —> D*^K±) = ( 16.3 ± 1.2 (stat) Iqj (syst) ± 4.8 (norm) ) x 10 -5 , 

where the third uncertainty is due to the uncertainty on the branching fraction of 
the normalisation channel. 
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1 Introduction 


The weak phase 7 is one of the least well-determined CKM parameters. It can be measured 
using time-independent decajj^] rates, such as those of B + D°K + or by time-dependent 
studies of B® — » D l f' >T K ± decays |1|. In time-dependent measurements with the decays 
h + , where h indicates a light meson, the sensitivity to 7 is a consequence 
of the interference between the amplitudes of the b —>■ u and b —> c transitions occuring 
through B® s yB^ mixing. The relevant Feynman diagrams for the B ° system are shown 
in Fig. [l] 



Figure 1: Feynman diagrams of the processes under study. The upper diagrams represent the 
two tree topologies (6 —> c and b —> u transitions, respectively) by which a B® meson decays 
into the D*^ 1 final state; the lower diagrams show the tree diagram of B® —> D*~tt + and the 
IT-exchange topology of B® —> D*~K + . 

The B° s —> DJK ± decay mode has already been used by LHCb to determine 7 with 
a statistical precision of about 30° |2|, in an analysis based on data corresponding to an 
integrated luminosity of 1 fb 1 . An attractive feature of B Q S —> D* S T K ± decays is that the 
theoretical formalism that relates the measured CP asymmetries to 7 is the same as for 
B® —» Df K ± decays, when the angular momentum of the final state is taken into account 
in the time evolution of the B®-B® decay asymmetries. 

The observables of the decay B® —» Di*'^K ± can be related to those of B° —> D^~ 7 r + 
as described in Ref. |Tj through the U-spin symmetry of strong interactions. This opens 
the possibility of a combined extraction of 7 . In addition, there is a higher sensitivity to 

1 Charge-conjugate states are implied throughout. 
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7 in Bg —* K ± decays than in B° —* 7%*) '7r + decays due to the larger interference 

between the b —> u and b —> c amplitudes in the former. 

The ratio 1Z = B{B—> DJK^/B^B® —» D~n + ) has recently been measured by 
LHCb [31 to be 1Z = 0.0762 ± 0.0015 ± 0.0020, where the first uncertainty is statistical and 
the second systematic. This is compatible with the predicted value of 1Z = 0.086±q;qq 7 from 
Ref. 111, which is based on SU(3) flavour symmetry and measurements from B factories. 
Under the same theoretical assumptions, the ratio 7 Z* = B(B ° — y D* T K ± )/B(Bg — y 
D*~tt + ) is predicted to be 1Z* = 0 . 099 i£o 36 Jl| and it is therefore interesting to test this 
prediction for vector decays. 

The B® — y D*~ir + and B® — > D* T K ± decays are experimentally challenging for detectors 
operating at hadron colliders because they require the reconstruction of a soft photon 
in the D*~ —>■ Dj'y decay. This paper describes the reconstruction of the B® —> D*~n + 
decay, previously observed by Belle (4] , as well as the first observation of the B Q S — * D* T K ± 
decay and the measurement of 7 Z*. This is the first step towards a measurement of the 
time-dependent CP asymmetry in these decays. 

The pp collision data used in this analysis correspond to an integrated luminosity of 
3.0 fb -1 , of which 1.0 ftU 1 were collected by LHCb in 2011 at a centre-of-mass energy of 
y/s = 7TeV, and the remaining 2.0 ftU 1 in 2012 at i/s = 8TeV. 

The ratio of branching fractions for the decays B® —> D* T K ± to B° s —> D*~n + is 
evaluated according to 


N K ± £ n + 
N n + £ K ± ’ 


( 1 ) 


where £x and Nx are the overall reconstruction efficiency and the observed yield, respec¬ 
tively, of the decay mode, and A" represents either a kaon or a pion (the “bachelor” hadron) 
that accompanies the D*~ in the final state. 


2 LHCb detector 

The LHCb detector [5J 6 ] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < p < 5, designed for the study of particles containing b or c quarks. The 
detector includes a high-precision tracking system consisting of a silicon-strip vertex 
detector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations 
of silicon-strip detectors and straw drift tubes placed downstream of the magnet. The 
tracking system provides a measurement of momentum, p, of charged particles with a 
relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200GeV/c. The 
minimum distance of a track to a primary vertex, the impact parameter, is measured 
with a resolution of (15 + 29/px) pm, where p x is the component of the momentum 
transverse to the beam, in GeV/c. Different types of charged hadrons are distinguished 
using information from two ring-imaging Cherenkov detectors. Photons, electrons and 
hadrons are identified by a calorimeter system consisting of scintillating-pad and preshower 
detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified 
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by a system composed of alternating layers of iron and multiwire proportional chambers. 

The online event selection is performed by a trigger which consists of a hardware stage, 
based on information from the calorimeter and muon systems, followed by a software 
stage, which applies a full event reconstruction. At the hardware trigger stage, events are 
required to have a muon with high px or a hadron, photon or electron with high transverse 
energy in the calorimeters. For hadrons, the transverse energy threshold is 3.5 GeV. The 
software trigger requires a two-, three- or four-track secondary vertex with a significant 
displacement from the primary pp interaction vertices (PVs). At least one charged particle 
must have a transverse momentum px > 1-7 GeV/c and be inconsistent with originating 
from a PV. A multivariate algorithm [7] is used for the identification of secondary vertices 
consistent with the decay of a b hadron. The px of the photon from D*~ decay is too low 
to contribute to the trigger decision. 

In the simulation, pp collisions are generated using Pythia |8j with a specific LHCb 
configuration (9|. Decays of hadronic particles are described by EvtGen 10 , in which final- 
state radiation is generated using Photos 
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The interaction of the generated particles 

with the detector, and its response, are implemented using the Geant4 toolkit |12 
described in Ref. 


as 
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3 Event selection 


Candidate B® mesons are reconstructed by combining a D*~ candidate with an additional 
pion or kaon of opposite charge. The preselection and selection for the two decays analysed 
for the measurement of 1Z* differ only by the particle identification (PID) [l4| requirements 
imposed on the bachelor tracks. The D*~ and D~ candidates are reconstructed in the 
D ~7 and K~K + tt~ decay modes, respectively. Each of the three D~ daughters tracks is 
required to have a good track quality, momentum p > lOOOMeV/c, transverse momentum 
Px > lOOMeV/c and a large impact parameter with respect to any PV. More stringent- 
requirements are imposed for bachelor tracks, namely p > 5000 MeV/c and px > 500 MeV/c. 
A good quality secondary vertex is required for the resulting D'-bachclor combination. 
Photons are identified using energy deposits in the electromagnetic calorimeter that are 
not associated with any track in the tracking system. Due to the small difference between 
the masses of the D*~ and D~ mesons, called A M in the following, the photons from 
the D*~ decay have an average transverse energy of a few hundred MeV/c 2 . A cut on 
a photon confidence level variable is used to suppress background events from hadrons, 
electrons and n 0 decays j 6 j. This confidence level variable takes into account the expected 
absence of matching between the calorimeter cluster and any track, the energy recorded in 
the preshower detector and the topology of the energy deposit in the electromagnetic and 
hadronic calorimeters. 

Additional preselection requirements are applied to cope with a large background mainly 
due to genuine photons that are not D*~ decay products, or hadrons that are misidentified 
as photons. The reconstructed mass of the D~ candidate and the reconstructed A m 


value are required to be in a ± 20 MeV/c 2 window around t-hcir known values 15 . The 
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Figure 2: (left) The K~K + tt~ invariant mass 
D*~tt + candidates. The points represent data, 
signal expected from the simulations. 



100 150 200 250 

A m [MeV/c 2 ] 


and (right) mass difference A m of the B® —> 
On the right plot the solid line represents the 


B», ^D-K+ (7r + ) decays are vetoed by a cut on the invariant mass of the D s K + 
(7r + ) system. PID requirements are applied to all final-state hadrons. Finally, the 
maximum distance in the rj-cp plane between the D~ and the photon is required to satisfy 
a JAt] 2 + A <p 2 < 1, where Ar) (A<p) is the pseudo-rapidity (azimuthal angle) distance 
between the corresponding candidates. 

To further reduce the combinatorial background while preserving a high signal efficiency, 
a multivariate approach is used. This follows closely the selection based on a boosted 
decision tree (BDT) 116,17 used in the measurement of the ratio of B® —> DfK ± to 


B^ D~^ 


branching fractions |3|. The algorithm is trained with simulated B° s —>■ D 


7T 


events as signal, and candidates in data with an invariant mass greater than 5500 MeV/c 2 
as background. The five variables with the highest discriminating power are found to be 
the Bg transverse flight distance, the photon transverse momentum, the Xip of the B° s 
candidate (where Xip is defined as the difference in y 2 of the associated PV, reconstructed 
with and without the considered particle), the angle between the B® momentum vector 
and the vector connecting its production and decay vertices, and the transverse momentum 
of the bachelor particle. Eight additional variables, among them the tra n sverse momenta 
of the remaining final-state particles, are also used. The trained algorithm is then applied 
to both the B° s —)■ D* s ^K ± and B° s —)■ D*~n + decays. 

The M(K~K + n~) and Am invariant mass distributions, as obtained from the decay 
mode B° s — > D*~n + , are shown in Fig. [2j These distributions have been obtained with all 
of the analysis requirements applied except that on the plotted variable. In both cases 
the B ° invariant mass is restricted to a ± 70 MeV/c 2 region around the known mass. A 
prominent peaking structure is observed in the Am distribution around 145 MeV/c 2 , due 
to the radiative D*~ to D~ decay. 
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4 Signal yields 


The signal yields are obtained using unbinned maximum likelihood fits to the B° s can¬ 


didate invariant mass distributions and are performed separately for B° s —y D* s 
B° s D*TK± decays. 


' 7 r + and 


The signal shapes are parametrised by a double-sided Crystal Ball (CB) function 18 


which consists of a central Gaussian part, with mean and width as parameters, and 
power-law tails on both lower and upper sides, to account for energy loss due to final-state 
radiation and detector resolution effects. The two mean values are constrained to be 
equal. When fitting the D*~n + and D* T simulated mass distributions all parameters 
are floated. When fitting data, the power-law tails parameters are fixed to the result of 
the fit to the corresponding simulation. Furthermore, both widths of the CB are set to 
those obtained from the signal simulation, scaled by a variable parameter in the fit to 
allow for differences in the mass resolution between data and simulation. The common 
mean of the double-sided CB is allowed to vary. 

Three background categories are identified. Partially reconstructed background decays 
are due to B® decay modes that are similar to signal but with at least one additional 
photon, as for example in the case of the B ° —>- D* T p± decays with p ± —>■ 7 r° (—> 77 ) n ± . 
Fully reconstructed background events are due to B° decays to the same final states as 
the B ° signal, D*~tt + and D* T K ± . The B° s —> D*~n + decays gives rise to a peak in the 
B ° — > D* T K ± decay mode when the 7r + is misidentified as a K + , a cross feed contribution. 
The cross feed due to K ± to ^ misidentification is negligible. Finally, a combinatorial 
background, where a genuine D~ meson is combined with a random (or fake) photon and 
a random bachelor track, can also contribute. 

The number of partially and fully reconstructed background components is different 
for each of the two final states. The invariant mass shapes for these backgrounds are 
obtained from simulation and are represented in the fit as non-parametric probability 
density functions (PDFs). The yields of these background components are free parameters 
in the fit, with the exception of the D*~7i + , D~p + and D*~p + contributions in the D*^K ± 
fit. The size of the D*~ti + cross feed is calculated from the D*~n + yield and the 7 r to K 

contributions are determined in a 


and D* p 4 


misidentification probability. The D~p^ 
similar manner, summed and fixed in the fit. 

To model the combinatorial background a non-parametric PDF is used. This is obtained 
from the events of the A m sideband in the interval [185,205] MeV/c 2 , with all other cuts 
unchanged. 

The results of the fitting procedure applied to the two considered decay modes are 
shown in Fig. [ 3 J The fitted yields are 16 513 ± 227 and 1025 ± 71 for the B ° —> D*~tt + 
and B ° —> D* T K ± cases, respectively. When the y 2 test is applied to gauge the quality of 
the fits, the latter fit has a y 2 value of 88.5 for 100 bins and 7 free parameters, the quality 
of the former fit is equally good. 

One of the distinctive features of the present analysis is the reconstruction of the decay 
mode D*~ —s- D ~7 at a hadron collider. The background-subtracted 7] and px distributions 
of these photons have been obtained using the invariant mass fit results described above 
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Figure 3: Invariant mass distribution of (top) B® —> D* 7r + and (bottom) B® -> D* T K ± 
candidates with fit results superimposed. The fitted signal corresponding to the first observation 
of B° s ->■ DpK* is shown by the dotted line in the lower plot. 


and the sPlot |19| method. These measured distributions are compared to the predictions 
of the simulation in Fig. |4j It is noted that most of the measured photons are very soft, 
with the average pt well below 1 GeV/c. 
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Figure 4: Distributions of (left) q and (right) pt of the photons for the D*~ir + (blue) and 
D* T K T (magenta) decays. Data, background-subtracted using the sPlot method, are represented 
by points, and simulations by solid lines. 


5 Systematic uncertainties 


Potential systematic uncertainties on 7 Z* are those due to the background modelling and 
the analysis selections, including the BDT and the PID cuts. Their effects are shown in 
Table [l] as relative variations of the final result, with their sum in quadrature assigned as 
the overall systematic uncertainty. The order in which the systematic uncertainties are 
described in the following text corresponds to successive rows in Table [lj 

Combinatorial background modelling uncertainties are studied by varying the default 
A m range used for the combinatorial background determination, [185,205] MeV/c 2 , to 
[205,225] and [225,245] MeV/c 2 . An alternative modelling of this background, using a 
parametric shape obtained from the D~ mass sidebands, is also tested. Finally, the 
statistical uncertainty due to the number of events in the range [185,205] MeV/c 2 is 
evaluated using the bootstrap technique 120,211. The corresponding uncertainty is taken 
to be the largest spread among the four differents checks. 

The uncertainty due to the finite size of the simulated samples used to study the 
partially reconstructed backgrounds is studied using the bootstrap technique. 

The uncertainties due to the D*~tt + cross feed and the D~p + and D*~ p + contributions 
to the D* T K ± fit are estimated by varying their expected yields. For the D*~ti + cross feed 
the ±lcr variation is obtained using the D*~tt + fit results. In the D~p~* 


and D* p "* 


cases 


the branching ratio uncertainties and photon kinematic distributions are different from 
the D*~tt + ones so the uncertainty in the yields are large. These yields are conservatively 
varied by ±50 %. The observed differences in the final result are assigned as the systematic 
uncertainties associated with these sources. 

The systematic uncertainty associated with the BDT is studied by reweighting the 
simulation to improve the agreement with data |3|. 

The 7r and K PID efficiencies used for the bachelor track have been extracted from 
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a D* + —y D°ir + calibration sample and parametrized as a function of several kinematic 
quantities of these tracks. The uncertainties in this procedure, propagated to the final 
result, lead to the PID systematic uncertainty. 

The systematic uncertainty from the hardware trigger efficiency arises from differences 


in the pion and kaon trigger efficiencies which are not reproduced in the simulation 22 


The uncertainty is scaled with the fraction of events where a signal track was responsible 
for triggering. 


Table 1: Estimated systematic uncertainties on TZ*. 


source 

relative variation (%) 

combinatorial background 

+ 4.7 

- 2.2 

simulation sample size 

±1.4 

D*~tt + cross feed 

±0.8 

Di*^p + “cross feed” 

+0 

- 1.6 

BDT 

±0.5 

PID uncertainties 

±1.0 

hardware trigger 

±1.0 

total 

+ 5.2 

- 3.5 


6 Results 


The ratio of branching fractions, measured in this analysis for the first time, is 

TZ* = B(B° S -± Df^K ± )/B(Bg -± D*~ tt+) = 0.068 ± 0.005 (stat) (syst), 

where the overall systematic uncertainty is mainly due to the uncertainty on the com 
binatorial background estimate. The result for TZ* ro 

.. J° a -)• D*- 

PID efficiencies. This factor is determined to be 1.095 ± 0.016 and is dominated by the K 
to 7r PID efficiency ratio. 

The measured value of TZ* is consistent with the theoretical prediction of TZ* = 
0.099^0030 11 1 , within the very large uncertainty of the latter. The theory is found to 
provide a good description of the measurements for both 7 Z* and TZ |3|. Other theoretical 


differs from the uncorrected B® —y 
D* T K ± to B2 —> D*~n + events ratio by a factor depending on the simulation and the 


predictions of 7 Z* have been published in Refs. 23-27 


Combining the measured value of TZ* with the value of B(B® —> D* 7r + ) obtained by 
Belle [4] leads to 

B(B Q S —> D= ( 16.3 ± 1.2 (stat) (syst) ± 4.8 (norm) ) x 10 -5 , 
where the uncertainties are statistical, systematic and due to the uncertainty on B(B° S —> 


D* s -7T- 1 
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